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Abstract 

The Hirota equation is a higher order extension of the nonhnear Schrodinger 
equation by incorporating third order dispersion and one form of self steepening 
effect. New periodic waves for the discrete Hirota equation are given in terms of 
ehiptic functions. The continuum hmit converges to the analogous result for the 
continuous Hirota equation, while the long wave limit of these discrete periodic 
patterns reproduces the known result of the integrable Ablowitz-Ladik system. 
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16 Conclusion! 5 

1 Introduction 

The Hirota equation jTS] 

iut + Ci (j)U:rx + q\ufuj + ic2 [pu^xx + 3g|u|^Mx) =0, i"^ = -1, C2pq ^ 0, (1) 

in which p,q,ci,C2 are real constants and m is a complex function of x and t, is quite 
important in the propagation of short pulses in optical fibers and lattice dynamics P^ 
123 1^. Indeed, this is one of the only two integrable cases |21l I2S] of a class of equations 

mm 

iut + aiu^^ + a2\u\^u + i [piu^^xx + (32\u\^Ux + /?3 (|m|^) uj = 0, (oij^Pj) real, (2) 

which incorporates the higher-order correction terms needed when the nonlinear Schrodinger 
equation (NLS) becomes insufficient or inadequate to describe the relevant physical phe- 
nomena. 

Just like the NLS, equation (P) is invariant under a Galilean transformation [2S1 QHl 
which could allow one to set Ci = 0, C2 = 1 without loss of generality. It is however 
preferable to retain both parameters Ci, C2 in order to better display the two limits, 

(a) the nonlinear Schrodinger equation (ci = 1,C2 = 0), 

iut+pur,:r + q\ufu = 0, (3) 

(b) the modified Korteweg-de Vries equation (mKdV) (ci = 0, C2 = 1,m = u), 

ut + puxxx + Squ'^u^ = 0. (4) 

If the time t is kept continuous and only the space variable x is made discrete, one can 
build a discretization of (^ which admits a discrete Lax pair [221 ^^^ therefore remains 
integrable. With the notation [/„ = U{X,t), X = nh, this semi-discrete equation is 

^Un,t + 2 [I'iUn+l + Un^l) + i/i(f/n+l " Un-l)] (1 + 6UjJn) ~ uU^ = 0, (5) 

and its continuum limit to (0) involves a convenient linear change of coordinates aimed 
at avoiding a Un^x term in the definition of (0), 

Un = U{X, t) = u{x, t), X = nh = X + Cot, 

6 = — h"^, u = 2cip/i~^, u = Qc2ph~^, cq = 6c2ph~^. ^ ' 

2p 

Bright and dark solitary (or localized) waves will occur when the signs of the disper- 
sive and nonlinear terms in the NLS portion of (jSj) are of the same and opposite signs 
respectively. The explicit expressions have been given earlier [19, 20 . In particular, bright 
and rational solitons on a continuous wave background are also obtained for the case of 
dispersive and nonlinear terms of the same sign. 

The main goal of the present paper is to deduce new, special periodic solutions of ^. 
First we shall start with solutions of the continuous Hirota equation (0), some of which, 
surprisingly, have not yet been documented. The discrete analogue of these solutions will 
yield exact results for (0). 



2 Continuum and discrete one-soliton solutions 

We begin by making a rather puzzling remark on localized pulses of these higher order 
envelope equations. In the defocusing regime {pq < 0), where conventional wisdom will 
dictate the presence of dark solitons, a localized pulse which rises above the mean level 
of a continuous background can exist. This solution, 



u 



^iiKx-ut)^^ 




k sinh ka 



cosh k{x — ct) + cosh ka 



+ iY 



2p 



P P coth."^ ka\ o 

c = C2P [k'^ — 7:{k coth ka)'^ + 3K^ — 6Y^ j , arbitrary 



(7) 



{a,k,Y). 



can be obtained with the method in [71 App. A]. A discrete equivalent of (|7|) is not yet 
known. 



3 Elliptic travelling wave solution in the continuum 
case 

Doubly periodic solutions of the continuous Hirota equation ((H) have been documented 
earlier jHlini, they are analytically identical to those |S.. of the Sasa-Satsuma case [23j of 
the integrable higher NLS (i.e. 3/3ia2 = /52«i,/?2 = 2/?3 in Q). Most of these solutions 
can be expressed as either a purely real elliptic function, mulitplied by the sinusoidal (or 
exponential) phase of wave oscillations. We propose a new form of solution expressed with 
nonzero real and imaginary parts. We shall present this solution first in the notations 
of Halphen (Appendix A). Although this notation appears to be more abstract and less 
common, it enjoys the advantage of displaying clearly the ternary symmetry that exists 
in the complex plane between e.g. the three Jacobi functions (en, dn, sn). We then employ 
the more conventional Jacobi elliptic functions. 

Assuming a solution of (0) in the form of a traveling wave (see Appendix A), 



u = ai 

one finds the unique solution 
2__2p 



haix - ct) +ibi 



hpjx - ct) 
h^{x - ct) 



^i{K{x—ct)—u)t) 



uj/p = ci{2hi - 2biK ~K' - 3e„) + C2(2fr(3e<, + K') + 661 (e^ - e^) + QK'), (9) 
c/p = 2ci{bi + K) + 3c2(e„ - K^ - 2biK - 2bl), 



. bi{e^ 



&i)(e^-e/3) = 0, 



which, as displayed in the last line, splits into three subcases 
(61 = 0): the same elliptic solution [H] as for the pure NLS case, 

{e-y — Ca — bl = 0): another elliptic solution, different from the previous one. 



{e.y — Cfs) = 0: a trigonometric solution, identical to the well known dark one-soliton 



u = ao 



— tanh — (x — ct) + ibi 



In each case, the solution depends on three arbitrary constants, e.g. respectively, (e^, Cq, K), 
(cq, bi, K), (cq, bi, K), i.e. three less arbitrary constants than the general travelling wave 
solution whose expression is still unknown. This lack of arbitrary constants is responsible 
for the factorized form of the last line in Q . 

To the best of our knowledge, the elliptic solution with 61 7^ (second case) is new. 

In the more familiar Jacobi notation [17 , this elliptic solution for 61 7^ is a periodic 
oscillation of amplitude Ai, 



u = Ai 



cn(r(a; — ct), k) + iB 



sn(r(a; — ct), k) 
dn(r(a; — ct), k) 



^i(K(x~ct)—ujt) 



;ii) 



where the wave vector r, the Jacobi modulus k, the velocity c and the angular frequency 
a; are related by 



f P = 1 - Bl 

Al = %\l-Bl), 

c/p = 2ci{K - rBi) + C2(r2 - 3K^ + GKrBi - Sr^B^), '^^^^ 

00 /p = ci{-K^ + 2KrBi + Ar^B^ - r^) 

+ C2i2Kr^ + Qr^Bi + 2K^ - QK^rBi - AKr^B^), 

and the three parameters {ax,Bi,K) are arbitrary. 

Remark. The nonintegrable higher NLS equation Q also admits the solution (jHJ, 
however with one less arbitray constant since K is then fixed. 



4 Particular travelling waves in the discrete case 

If one assumes for (0) the same kind of solution as for (^ , namely 



f/ = ai 



ha{X -ct) +ib_ 



hfsjX - ct) 
h^{X-ct) 



^i{K{X-ct)-u)t) 



(13) 



one finds solutions similar to those of the continuous case. In particular, the elliptic 
solution with 61 7^ is 



2__h^(M _ 2 



n 



5N^ 



uj = V \l 7^^ + /iOi(e/3 - e^j 



iVa 



A^. 



(14) 



^ ^ ^^^^^ /ih.(/.)h,(^) + z.&,h,(/.) ^ ^^^.^^^^^ ^ ^^^^ ^^^ ^^^ 

i Vrv 



in which N^ denotes the nonzero constant 

iVa = \ii{h) - (e„ - e/3)(e„ - e^) = 2ha(/i) \ip{h) h^(/i) ha(2/i). (15) 

The corresponding expression in the Jacobi notation is 

sn(r(X — ct), k) 



U = Ai 



cn{r{X - ct),k) + iBi 



dn(r(X-ct),fc) 



^i(K(x—ct)—ujt) 



(16) 



in which the relations linking {r,k,c,uj, K, Bi,Ai) are easily deduced from ()14|) using the 
correspondence 



5 Special limits 

5.1 The continuum limit 

The continuum limit h ^ is easier to perform on ^^. The fixed constants (5, yU, z/) 
behave as indicated in ©, the elhptic functions of h expand as 

h„(/i) = h-' - —h^ + OQi"), N^ = h-^ - 2e„/i-2 + 0(/i°), 



iV„ \ 2 2 

h«(/i)h^(/i) _ ^ ^^ 3e^ + 5e^ ^^ 



hii- " ^ h' + 0{h 



therefore 



c = Qc2ph-^h (l - ^^^±^/,2 + o{h^)\ + 2c,ph-\h\ (18) 

X — ct = a; — (c — co)t = a; — p(2ci6i — 3c2(3e/3 + 5e^))t + 0{h^), 
and the solution (fT^ indeed goes to the solution Q. 

5.2 The long wave limit 

In the long wave limit k -^ 1, the solution (fTH|) goes to JH! 

U = a2^^ sech k{X - ct)e^^^''-'\ 4 = ^, 
h q 

uj = u + {fi sin Kh — v cos KK) cosh fc/i, (19) 

, ^^, . ^^, , sinh A;/i . 

c = [fi cos Kn + usmKhj ^ , arbitrary = [k,K). 

K 

6 Conclusion 

Periodic waves of the discrete Hirota equation have been generated from their counterparts 
for the continuous Hirota equation. 

Equation © combines the main features of the integrable discretization of NLS |21 El 

ng, 

^-^ + J^2 + (^"+1 + ^n-lli'n^n = 0, (20) 



and the discrete modified Korteweg-de Vries equation [T^ 

^ ^ (V'n+l - V^n-l)(l+^n) 
9t ~ /l2 



(21) 



The most general travelling wave solution of the continuous Hirota equation is an 
hyperelliptic function of genus two, as detailed in Appendix B. The elliptic solution 
is a very particular case of this general traveling wave. Future work should address the 
question of incorporating the freedom of the three missing arbitrary constants in the 
present elliptic solution. 
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Appendix A. The basic elliptic functions 

In order to fully display the ternary symmetry of the Jacob! notation, and consequently 
make the practical computations much easier, following Halphen ^T], we will choose the 
basis functions as 



ha{u) = ypiu) - e„, a = 1, 2, 3, (22) 

in which p is the Weierstrass elliptic function defined by 

p'^ = 4(p - ei){p - e2){p - es) = Ap^ - g2p - 93- (23) 

The link with more usual (but less symmetric) notation such as (cs, ds, ns) or (sn, en, dn) 
is classical H §16.10, 16.20, 18.9.11], e.g. [IH Chap II, (16) p. 46], 

cs(^) ^ ds{z) ^ ns{z) ^u^ 1 
hi{u) h2{u) hsiu) z y/ei - 63 ' 

and our main results will also be displayed in this more usual notation, familiar to the 
physicists. 

Since the three factors p — ej,j = 1,2,3 of p' are indistinguishable, any solution 
depending on (hi,h2,h3) remains a solution after any permutation (a,/5,7) of (1,2,3). 
Therefore, in the complex plane, a hi solution represents any of the twelve familiar func- 
tions pq of Jacobi (p,q among s,c,d,n). However, on the real axis of x or of x — ct, 
depending on the values of (ci, 62, 63), less than twelve functions will be admissible. 

For deriving the continuum limits, it is useful to consider the Laurent expansion at 
the origin [HI Chap VII p. 237] 



h<^(«) = - - ^" 2 + y'' - ^'-) 40 + 2^.5.7 
, 28^1 - 225e«^3 - 105e2^2 7 



u' 



29.3.52.7 



u' + 0{e). (25) 



Full details on this symmetric notation of Halphen can be found in Ref. [Hj. 
In the main text, (a, /5, 7) denotes any permutation of (1, 2, 3). 



Appendix B 

The traveling waves of (P) are defined by tlie reduction 

in wfiicfi M , ^p are real functions and c, uo real constants. 

The two real functions (M, Lp) obey the sixth order differential system 



(26) 



C2 



+iC2 



+Cl 



3(/?'M" 3(^'M'^ 3(^"M' 



2M 



+ 



/// , /3 3g<^'M 



4M2 2M ■ p 

3v9'^M' 3gM' 



2M 4M2 



3M'M" 3M'^ 
+ 



2M AM- ^ +p^ 



8M3 2M 



^ 2p 



- 3^V 



cu)' + UJ . cM' 

H 1 

2pM 



0. 



P 



(27) 



All the first integrals of the system (j7fj) can be generated systematically from the Lax 
pair pP with the result, 



K, = cl 



M" M'^ 



/2; 



+ TT7 + f M 



AM 



qM"' 



2 , , . cM 

-CiC2^M + C2-— , 

3 3p 



(2^ 



^' /\^'2 



A2 = C2 -V? iW + 



2M 



+ p:"M' + 2(^'"M 



K. 



+ CiC2(M" 



M'^ gM2 



-'2; 



AM 



M 
M'^M" 
AM- 



P 



Aifi' M + 2ciip'M - (cic + C2Uj[ 



M 
P 



(29) 



_^'-M" + ^MM"+ ^' 



/4 



P 



16M3 



^ + 2ip'(p"M' + if' M + ^ 



-M^' 



2M 



2p 



pz 



^//\^'2 



+ C1C2 lif'M" - ^^^ - y^"M' - 2y;''M 



2M 



- ^M" + 
P 



C\ C2C 



A^ 2p 



' ^% A^'^ m] - 2''' + ''''-' 



M 



P 



if'M 



+ 



cfp — AC2C 



qM^ + 



CipUJ + c 
p2 



M. 



(30) 



The associated spectral curve, derived from the monodromy matrix M. 
dei{M- li) = li- -P{\) =0, 
P(A) = - {ac2p\^ + 2cip\^ + c\ 



f)-pq[QK,X--K2X+'^ 



(31) 



is hyperelliptic and it has genus two. Therefore the general solution is a singlevalued 
function of ^ expressed with genus two hyperelliptic functions. To our knowledge, neither 
the separating variables nor the explicit expression of this general solution are known. 
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When the polynomial -P(A) has a double root, the solution is elliptic and our goal here 
is to derive explicitly this elliptic travelling wave solution. It depends on the following ar- 
bitrary constants: the origins ^o, V^o of ^ and ^9, plus four constants out of c, u;, Ki, K2, K3. 
After this is achieved, we want to derive the corresponding solution of the discrete equa- 
tion (jnj just by the natural discretization x — > n/i of this elliptic solution of the continuum 
equation. 
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